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We discuss the relevance of the improved mean-field slave-fermion theory to describe the Nagaoka 
{U = oo) limit of the Hubbard model. In this theory the crucial on-site constraint of no double 
electron occupancy is taken into account rigorously prior to the mean-field approximation. At one- 
loop approximation the effective mean-field action shows a long-range ferromagnetic order over the 
whole doping range. This indicates that the slave-fermion mean-field theory does not constitute an 
appropriate framework to describe the physics of the Nagaoka phase. We discuss the drawbacks 
of this mean-field theory and present some results on the derivation of a low-energy effective spin 
action to describe the Nagaoka phase beyond the mean-field approximation. 



I. INTRODUCTION 

It is now widely believed that the essential physics of 
strongly correlated lattice electrons is encoded in the lo- 
cal constraint of no double electron occupancy. This con- 
straint prohibits double electron occupancy of a lattice 
site due to the large local Coulomb repulsion between 
hopping electrons. The on-site Hilbert space is thus re- 
stricted to states with at most one electron per site. Such 
a modification of the underlying Hilbert space results in 
dramatic consequences for the low-energy properties of 
the relevant electron systems revealing a rich and un- 
usual physics in this limit. 

Formally, given a local electron operator with the 
spin projection tr =t,i, the local no double occupancy 
(NDO) constraint reads 



< 1. 



(1) 



Let us denote by Cia the projected electron operators 
that satisfy this requirement. Since the inequality ([T]) 
is hard to deal with analytically, one can try to circum- 
vent this difficulty by turning Eq.lT]) into an equality at 
the expense of the introduction of some redundant de- 
grees of freedom. The two well-known ways to accomplish 
this are the so called slave-fermion (SF) and slave-boson 
(SB) representations of the constraint electron operator. 
Namely, by introducing the "slave boson"—, one decom- 
poses the on-site constraint electron operator in the form. 



(2) 



where hi is a charged spinless boson (holon) , while fi^ ia 
a neutral, spin 1/2 fermion operator (spinon) satisfying 
the NDO constraint 



(3) 



Alternatively, one can also introduce a spinless fermion 
fi to describe the charge degree of freedom and a spinful 
boson hirj to keep track of the spin degree of freedom. 
This is the "slave fermion" approach-. 



The NDO constraint now reads 



1. 



(4) 



(5) 



The electron spin operator, Qi = 5 Z^o-o-' ^L^fo-'Cjo-' , 
takes in this representation the form 



(6) 



where Eq.([5|) has been used. Here T*'s are the Pauli 
matrices. 

There is an apparent U(l) gauge redundancy in the de- 
compositions of the projected electron operator as given 
by Eqs. (|2l4p . Namely, the local gauge transformation 
generated by the NDO constraint 



(7) 



keeps the representations ([2]) and ([4]) intact. The slave- 
particle gauge theory is infinitely strong coupled since 
there is no kinetic energy for the bare gauge field. Con- 
sequently, one should exercise some care in applying an 
ad hoc mean-field (MF) approximation to treat the slave- 
particle theories. Such an approach can only be justified 
provided that the low-energy gauge coupling gets effec- 
tively renormalized to some finite and, potentially, weak 
coupling. However, this seems hardly to be the case in the 
slave-particle theoriesS., which in turn casts some doubts 
on the results obtained within the MF slave-particle theo- 
ries. Indeed, both the SF and SB theories should in prin- 
ciple produce physically identical results for the t-J model 
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of strongly correlated electrons, however they give in the 
MF approximation very different phase diagrams^. 

The strong electron correlations are at work to full 
extent in the so called Nagaoka {U = oo) limit of the 
Hubbard model. Indeed, in this case an infinitely strong 
Coulomb repulsion strictly prohibits the double electron 
occupancy of the lattice sites, and the NDO constraint 
becomes of the utmost importance. In the infinite U 
limit, the Hubbard Hamiltonian reduces to 

^ = -Yl ^V^aCja + /i ^(1 - SLq.), (8) 

where tij is a symmetric matrix whose elements represent 
the hopping amplitude t between the nearest-neighbor 
sites and are zero otherwise. We have also introduced 
the chemical potential /i to control the total number of 
vacancies (holes). Despite its seemingly simple form, this 
Hamiltonian cannot be diagonalized due to the fact that 
the projected electron operators fulfill complicated com- 
mutation relations resulting from the explicit manifesta- 
tion of strong correlations. 

Since the Coulomb repulsion whatever large it may 
be, is in practice always finite, the one band U = oo Hub- 
bard model is basically a toy model and does not describe 
any specific material^. Nevertheless, it captures an ex- 
treme limit for the physics of strong electron interactions. 
Had the projected electron operators in ([5]) been replaced 
by their related conventional operators, the model would 
have been reduced to a system of nonintcracting elec- 
trons revealing thereby a trivial physics. On the other 
hand, the physics behind the model ^ is certainly far 
from trivial. Indeed, Nagaoka^ proved a theorem stating 
that for one hole the ground state of the U = oo Hubbard 
model is a fully saturated ferromagnet. This provides an 
interesting example of a quantum system in which ferro- 
magnetism appears as a purely kinetic energy effect with 
hole hoping (itinerant ferromagnetism) emerging as a re- 
sult of the strong correlations from the NDO constraint. 

Unfortunately, despite extensive work over many years, 
this model and itinerant ferromagnetism are still poorly 
understood. One of the important questions to be ad- 
dressed concerns the thermodynamic stability of the Na- 
gaoka phase. That is, whether or not the Nagaoka state is 
stable when the density of holes is finite in the thermody- 
namic limit. There are arguments both for— li^iiiii^ii^ii^ 
and againslfi^ii^iiiii^'i^ the thermodynamic stability of 
the Nagaoka phase and comparisons between various ap- 
proaches have been made carefully (for a recent example, 
see^^) . The basic problem that prevents one from reach- 
ing a definite conclusion on that is the large- [/ limit or, 
equivalently, the local NDO constraint which is very diffi- 
cult to deal with in a controlled way. For example, the nu- 
merical studies of this problem, although well developed 
for the Hubbard model, cannot incorporate the large-C/ 
limiij^. Analytical approaches, on the other hand, ba- 
sically imply a MF treatment in which the local NDO 
constraint is uncontrollably replaced by a global condi- 
tion. 



In this sense, it is interesting to note, that 
there are several works considering different variational 
approache a^i^^'^^'^-^1^^'^'^1^^ in which the NDO constraint 
is automatically built into the trial state vector. In gen- 
eral, the following observation holds: variational esti- 
mations that involve more realistic refined trial wave- 
functions result in a smaller value of the critical hole 
concentration. 

One may therefore think that a straightforward appli- 
cation of the MF theory that treats electron correlations 
on average is not likely to produce reliable results. Specif- 
ically, in the present paper we show that the treatment 
of the Nagaoka phase within the MF SF approximation 
results in the qualitatively incorrect conclusions. 

We outline below the main motivation of this work. 
The standard SF MF theory which treats the NDO con- 
straint only on average is known to predict a stable FM 
phase for the U — oo Hubbard model over a finite though 
unphysically large doping range. Since the local electron 
correlations encoded in the NDO constraint are surely 
of the utmost importance, if one attempts to deal with 
the properties of this model in a correct and reliable way, 
one would naturally think that a proper treatment of the 
constraint prior to the MF approximation could improve 
that MF result shifting the critical hole concentration 
towards a much smaller and physically reasonable value. 
In this paper we show that this is not the case. On the 
contrary, the improved MF theory predicts a stable FM 
phase over the whole doping range. 

Two main routes which redirects of this result are dis- 
cussed in Sections II and HI. We argue that the finite 
critical hole concentration predicted by the standard MF 
theory is just an artifact of the uncontrolled treatment of 
the NDO constraint. Moreover, the theory based on the 
MF treatment of the spin degrees of freedom affects the 
physics of the problem in a qualitative way. Not only the 
NDO constraint needs to be fulfilled. It is clearly very im- 
portant to treat the electron correlations encoded by that 
constraint in a proper way considering both charge and 
spin degrees of freedom on equal footing as dynamical 
rather than the MF variables. Some preliminary results 
in this direction are reported in Section IV. 

II. MF SF TREATMENT OF THE NDO 
CONSTRAINT 

In this section we briefly go over the conventional MF 
SF theory of the ferromagnetic phase in the U = oo Hub- 
bard model as given by, e.g., Boies et a&^. 

To set up the stage, let us start with the SF repre- 
sentation of the projected electron operators given by 
Eas. (|4l5[ ). The Hamiltonian ^ can be then rewritten in 
the form 

H = Y.^^f]blb,.f^+^^Y.fU^^ (9) 

where Eq. (O has been used. The partition function 
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takes the form 

Z = f DXD4>D'il;Y[DzaDzae^^''''''^\ (10) 

with the action 5* = J^^ dTCsFij) and the Lagrangian 

^Sf{t) = ijjj [{-dr - /i + - U] Z^„Zj„]%l)i 

ij a 

i(7 i 

Here z^g. and ip^ stand for complex numbers and complex 
Grassmann parameters, respectively. The purely imag- 
inary field Xi has been introduced to enforce the local 
NDO constraint, 

V'iV'i + X! ^i'^^i'^ = 1- (12) 
a 

In this representation the electron spin operator Q be- 
comes 

A. Zeroth-order approximation 

At the zeroth-order MF approximation boson variables 
take a macroscopic value and all fluctuations in space and 
imaginary time are neglected, i.e. 

z,.(t)=40), A,(t) = A(o). (14) 

Under this condition, the fermionic degrees of freedom 
in can easily be integrated out to yield for the free 
energy := -^logZ, 

= -]^Eiog(i + ^"'^'^^''^""'"'^) 

- ^|40)pA(0)+A(")-M^, (15) 

a 

where 

i?(fc) = (^|z(")|%^ 

and (5 is the average density of holes. The ground-state 
saddle point equations determining li and A^o) de- 
scribe a saturated ferromagnet, m oc (1— (5), where the av- 
erage magnetization is given by m = ^ crlz^^'' p. Such 
a solution exists over the whole doping range, < (5 < 1. 

It is interesting to note that due to the global MF 
constraint 

5 + 5:4°M°) = i, 



the spin variables and z| cannot vanish simultane- 
ously. This observation rules out the possibility of hav- 
ing, within that approach, a phase transition into a dis- 
ordered phase by identifying the set of macroscopic spin 
amplitudes, z^a^ ^ as the relevant order parameter. 



B. One-loop approximation 

The one-loop correction to the effective action can in 
principle alter this result destabilizing the zeroth-order 
MF theory at some critical hole concentration, 5c- In 
the MF SF theory the one-loop approximation takes the 
forir^: 

z,^ = z("'<5,,o + (5<, (16) 



Ag = A(o)(5,,o + '5A„ (17) 

where the zeroth-order MF solution, z^^^ — and 

z^^^ — \/\ — Equation (|16p tells us that quantum 
fluctuations of the magnetic "order parameter" are sup- 
posed to be small. Besides that, Eq. p7|) implies that 
the constraint-generated gauge field is slowly fluctuating 
function around its MF constant value. Assuming this, 
one can substitute Eqs. (|16ll7p into the SF action to 
perform the fermion integral (fTO|) at one-loop level. In 
this way, one arrives at a purely bosonic effective action. 
Analyzing then the stability criterion for the remaining 
bosonic integral over 8z„ and 5A, one can evaluate the 
critical value, be ~ 0.7, at which the zeroth-order MF 
theory becomes unstable against quantum fluctuations. 

However, this value of be is too much large to agree 
with the recent results of the variational Monte Carlo 
studies^^ ((5c = 0.38) and small cluster calculations^ 
((5c = 0.22). Note that the smallest value of (5c obtained 
so far is (5c — 0.17>i^. This MF value also contradicts the 
exact result obtained in infinite space dimension^^. 

The one-loop MF result can be improved, in princi- 
ple, by imposing the crucial NDO constraint prior to 
the MF calculations. Surprisingly this is not the case. 
On the contrary, as we show in the next section an im- 
proved MF one-loop theory predicts exactly the same 
result that follows directly from the zeroth-order MF ap- 
proximation: the FM phase is stable in the whole doping 
range, < (5 < 1. This apparently indicates that the MF 
theory based on the MF treatment of the spin degrees of 
freedom affects the physics of the problem in a qualita- 
tive way. Within that approach, the emergent long-range 
ferromagnetic order appears to be just an artifact of the 
MF approximation, rather than a dynamically generated 
effect. 
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III. SLAVE-FERMION EFFECTIVE ACTION 
WITH THE EXPLICITLY RESOLVED NDO 
CONSTRAINT 



where 



Let us turn back to the exact representation ((10)) . As a 
first step to proceed, we rigorously resolve the constraint 
V'iV'i + X^cr ^i'^-^io' ~ ^- a rcsult, we arrive at the 
SF path-integral representation of the partition function 
that explicitly incorporates the NDO constraint. 

The constraint can be explicitly resolved by making 
the identifications 



1 + ZiZi+ ^j^i 

l + ZiZi -I- CiC. 



1 + ZiZi+ ^j^j 



(18) 



with the variables 2^1,^1 and 1/)^ being free of any further 
constraints. 

Note that the local gauge transformation ^ reasserts 
itself in the form, (f>i — > 4>i + 9i. In contrast, the projected 
variables 

Zi^ZiJZil^, S_^^1pi/Zi^ 

are seen to be manifestly gauge invariant. As we al- 
ready mentioned, this gauge symmetry is a consequence 
of the redundancy in parametrizing the electron opera- 
tor in terms of the auxiliary boson/fermion fields. The 
gauge ambiguity related with the redundancy of the SF 
representation is now expressed by a single variable (pi. 

The domain of the flat measure in pUj) that involves 
the spin up bosonic fields can be rewritten at every lattice 
site as Dzi^Dzi-^ = D\zi-^\'^D(j)i. The jzi-f P field can easily 
be integrated out from Eq. fTU]) thanks to the constraint 
(fT2)) . Since the action in (fTO| is U(l) gauge invariant and 
hence independent of ipi, the integration over 4>i results 
merely in the appearance of some numerical factor (vol- 
ume of the gauge group) that can be taken care of by a 
proper normalization of the partition function. For the 
remaining integration we have (the site dependence for 
the moment being suppressed), 

Dz^Dz^Di^Di; = sdet|| ^^7'^i' ^'^^ || dzdzd^dl 

The Jacobian of this transformation (superdeterminant) 
was evaluated in^* to yield 



sdet 



| ^(ai,Qi,/,_/) | 
d{z,z,^,Q 



1 + + 



Putting everything together we get a new represen- 
tation for the SF partition function IjlOp with the local 
NDO constraint being explicitly resolved to give 



J Dt,{z,0 e 



(19) 



3,i 



2TTi l+kjf+CiCj 



stands for the measure with the boundary conditions, 
2:^(0) = Zj{j3), = ^^jiP)- The action now reads 



Jo 



H^'dt. 



(21) 



The first part of the action (PT|) is a purely kinematic 
term that reflects the geometry of the underlying phase 
space while the classical Hamiltonian becomes, 



^ 1 



-f ZjZi) -I- H.c. 



(l + kp+c.c.)(l + kip+^.-c,) 



(22) 



The new set of the gauge invariant variables, (2,^), ex- 
plicitly resolve the NDO constraint at the apparent ex- 
pense of a more complicated compact phase space for the 
projected electron operators. 

Finally, we make a change of variables to decompose 
the full measure given by (|^D|) into the product of the 
conventional spin and fermion measures, 

i^M.p»(^,^] - ll27r*(l-f|z,(t)P)2' 

respectively. Such a reparametrization can be taken to 
be 



(23) 



Up to an inessential factor which redefines the chemical 
potential, we get 

«,e), (24) 

and the effective action becomes 
S = ^ ia.i{T)dT - ^ / £i {dr + l-J. + ia^) £,idT 

i Jo . "'0 



/3 



H^^dT 



(25) 



This action involves the U(l)-valued connection one-form 
of the magnetic monopole bundle that can formally be 
interpreted as a spin "kinetic" term, 

,a^-{z\d,\z)^Y^—^, 
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with \z) being the su(2) coherent state (see Appendix). 
This term is also frequently referred to as the Berry con- 
nection. The dynamical part of the action takes the form 



(26) 



Here {zi\zj) stands for an inner product of the su(2) co- 
herent states, 



{z,\Zj) 



1 + ZiZj 



The classical image of the on-site electron spin opertor 
(fT3)) reduces to 

where 5^' is given by Eqs. ([5T|) . 

The action ([^5]) is invariant under the global SU(2) 
rotations that now take the form 



UZi + V 



e<->e^^'6, a,^a,-dC^, (27) 



-VZi + u 



where 

= -ilog 



'VZj + u 



-VZi + u 



e SU(2). (28) 



where 



(/)(t) = — — / ads + I ads. 
P Jo Jo 

Note that (^(0) = <j>{f3). The effective action then becomes 

/ ia{T)dT~ / ^{dr+fl)^dT, 

Jo Jo 

where M = + ^ /q^ iadr. The partition function is now 
given by the path integral 



where the measure factor is represented by Eq.(24). In- 
tegrating out fermions yields 



Z 



Since 



and J Dfispin — 1, one finally gets Z — 2 + e^''^, which 
is indeed what it should be desired. 



Equations (25,26) provide a rigorous representation of 
the effective action of the U — oo Hubbard model. One 
can therefore wonder if these equations lead to the ex- 
pected result if one nullifies the tij hoppings. In this 
limit the infinite U Hubbard model (8) reduces to the 
exactly solvable reduced Hamiltonian: 



i/ = A*^(l-cLc.,)=^X 



00 
i 5 



where the Hubbard operator X^^ is represented by the 
diagonal 3x3 matrix with eigenvalues 0,0 and 1. As a 
result, the partition function reduces simply to: 

Z ^2 + 3-1^^. 

As a check of the validity of our representation our equa- 
tions (|25l26p should also recover this result. 

To see that this is indeed the case, let us consider the 
on-site action (25) at < = 0, 



S 



P rP _ 

ia{T)dT —I ^ [dr + fi + ia) ^dr. 
Jo 



By a gauge transformation the potential a(T) can be 
brought into a time independent form, 
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A. Improved MF theory: zeroth-order 
approximation 

In this approximation the bosonic spin variable Zi(t) 
takes on a macroscopic value z*^"^. Since (z^'^^jz^^)) = 1, 
the spinless fermion MF Hamiltonian reduces to the 
representation 



(29) 



This zeroth-order Hamiltonian does not depend on z^^^ 
whose value merely determines a direction of the total 
electron magnetic moment. 



QMF = ^e/(z("),z(0))(l-^), 



(30) 



with the explicit representation of the spin moment 
Sci{z, z) being given in Appendix. 

^From now on we take the total electron magnetic 
moment aligned along the z-axis. To achieve this, one 
needs to set z(°) = (see Appendix). The Hamiltonian 
([29|) along with Eq. ([30|l describes a fully polarized fer- 
romagnet for any hole concentration 5 < 1. Therefore, 
the MF theory that takes care of the NDO constraint at 
the outset, produces naturally, at zeroth-order approx- 
imation, nearly the very same result that is produced 
by the zeroth-order MF theory (fH|) which only treats 
the NDO constraint globally. Note, however, that in the 
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representation ((29)l . in contrast to Eq. pS)) . there is no 
renormalization of the fermionic bandwidth. Again, one 
sees that due to the NDO constraint, the MF treatment 
of the spin dynamics, Zi{t) = z^^\ automatically drives 
the system into an ordered phase. Regardless of any spe- 
cific value of the constraint- free spin variable z^^\ the 
system always stays in the ordered ferromagnetic phase 
as dictated by Eq. ([50]) . The "order parameter" z^-^\ 
once again, fails to describe a phase transition out of the 
ferromagnetic phase. 



B. Improved MF theory: one-loop approximation 

In the present subsection we derive the one loop ap- 
proximation for the effective MF spin action, with the 
local NDO constraint built in from the outset. 

To this end we take into account the Gaussian fluc- 
tuatins of the bosonic spin variables around their MF 
value, 



(31) 



and expand the action (j25p up to the quadratic order in 
the new 6zi variables. The fermionic path-integral can 
then be evaluated to this order explicitly and we end up 
with the following one-loop spin MF effective action, 

S ^ So + AS, (32) 

where 5*0 represents the zeroth-order MF action, whereas 

AS* = TrSiia) + TrG(^t)ST.. (33) 

Here 

S{iai) = ^{SziSzi - SziSzi) 
is a linearized spin kinetic term, 

6T.ij = -Uj{SzjSzi - ^\Szi\'^ - i|(5zjp). 



and 



^{t) = (3t - M - k 



is a zeroth-order MF Green function. The trace has to 
be carried out in both space and time indices. 
In the momentum space the action p3p reads 



AS" = / {SzqSzq- 6z^6z^)dT 



where 



N 



(34) 



(35) 



and ff; = (^e^(*p+t^) + is the Fermi distribution func- 
tion. The action ([34|) corresponds to the bosonic spin- 
wave Hamiltonian, 



(36) 



with the ferromagnetic spin-wave dispersion relation, 
w^cx , g -> 0. 

These quantum fluctuations cannot destabilize the 
zeroth-oreder MF solution, since the excitation spectrum 
LUq' is a non-negative function of the hole concentration S, 
provided tp — t-p. This can be proven in the following 



way 



.29 



^ fp{^q+p V) 

peBZ 

U E /.-(w-v) + i4 E /p(w-v) 



27V ^ 

peBZ 



2N ^ 

peBZ 



(37) 



In the second summation we substitute p — s- —p— q. In 
this way we get 

'^9 ~ oT? E/ fpi^q+p~^p) 



2N 
1 1 



peBZ 



2 j\j E/ f-p-qi^q-p-q ^-p~q) 



peBZ 



2 jY E/ fp(^Q+p V) 
peBZ 

9"a7 E/ fp+gi^p ^ ^p+q) 



2N 
1 1 



peBZ 



9 AT E/ ^-^P fp+q)i^q+p V) 



27V 

peBZ 

Since is monotonically decreasing function of tp, the 
quantity (/p-— fp+q){tq+p ~ tp) is always non- negative. 

In fact, the spin- wave stability condition, uj^ > Q, q ^ 
0, is only a necessary condition for the saturation of 
ferromagnetisni^. However, it is possible to show that 
if the semiclassical analysis is assumed to be true, this 
also becomes a sufficient condition. Therefore, if the MF 
(semiclassical) decomposition of the spin amplitude ([3T|) 
is taken for grant, the ensuing SF MF theory is spin- wave 
stable over the whole doping range. 

There is, however, a rigorous result that states that 
for a large enough hole concentration the Nagaoka state 
posesses an instability^i. This result is heavily based 
on the consideration that the spin degrees of freedom in 
the U = oo limit are in fact dynamical variables rather 
than nearly frozen spins in the background. As was dis- 
cussed above, the MF treatment of the spin degrees of 
freedom automatically drives the system into the ordered 
FM phase. Although in this case the SF MF theory works 
well, it is applied to a substantially altered Hamiltonian 
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that does not seem to bear much in common with the 
original problem. 



IV. ONE-LOOP APPROXIMATION BEYOND 
MF THEORY 



In the preceeding section we derived the improved one- 
loop MF theory to treat the Nagaoka's phase. In the 
improved theory we imposed the NDO constraint rigor- 
ously from the beginning and, only after that, the MF 
approximation was applied to describe the dynamics of 
the spin degres of freedom. The improved theory predicts 
a stable FM phase over the whole doping range. At the 
same time, the standard SF MF theory that treats the 
NDO constraint globally at the zeroth-order and locally 
at the one-loop order tells us that the FM phase is stable 
over a finite though quite large doping range. Since the 
improved MF theory should in any event provide a bet- 
ter description, one may conclude that the predictions of 
the standard MF theory are not reliable and are just an 
artifact of the uncontrolled MF treatment. 

In this sense, it is interesting to note, that in an ear- 
lier paper— a variational principle was formulated in such 
a way that the NDO constraint was automatically built 
into the trial state vector which considers that all the 
electron spins except one are aligned. The remaining 
electron is kept frozen. It can then rigorously be shown 
that at f/ = CX3 the Nagaoka ferromagnet is always sta- 
ble for any 6. This qualitatively agrees with the result 
following from the improved one-loop MF theory. 

However, if the remaining electron is also allowed to 
hop around, the FM ground state immediately becomes 
unstable for large enough doping^^. All this indicates, 
that even the improved MF approximation is too restric- 
tive in the sense that it qualitatively affects the physics 
described by the model by unnecessarily freezing the spin 
degrees of freedom. This leads to the physically incor- 
rect conclusion of predicting the FM phase for the whole 
doping ranged. It is evident that not only the NDO 
constraint needs to be fulfilled. It is also very impor- 
tant to treat the electron correlations encoded by that 
constraint in a proper way considering both charge and 
spin degrees of freedom on equal footing as dynamical 
variables. It would therefore be very appealing to ad- 
dress the problem of the thermodynamic properties of 
the Nagaoka phase starting right from the low-energy ef- 
fective quantum spin Hamiltonian rather than from the 
quasi-classical MF spin effective action. Technically this 
approach is, however, quite complicated. In view of that 
we report in this section only some preliminary results. 

Specifically, we derive the contribution to the low- 
energy effective spin action of the U = oo Hubbard 
model up to the lowest nonvanishing order in the spin 
self-energy = —t{zj\zi) beyond the MF approxima- 



tion. To this end, we rewrite action (1251) in the form 



S = i^y^ / CLi{T)dT 

+ Ur)G-^^{T,<j)^A<^)dTda, (38) 

where 

G'~-^(T,cr) = G^„j^^.(r,cr)-mi(T)%(5(T-(T)+Sy(T)(5(r-cr), 
with Tiij = —tij{zj\zi) and 

The fermionic degrees of freedom can formally be in- 
tegrated out to yield 

j D^Di exp (t)G-.i (t, a)e, {a)dtds 

= exp Tr log G"^ 

= exp (rr logG(o5 + Tr\og{l - G(o)m + G(o)S));39) 

Here the trace has to be carried out over both space and 
time indices. Calculating explicitly the contribution com- 
ing from the zero order Green's function, it yields: 

Zo := Zt=a=o = cxp(rrlogG^p 



2cosh — 
2 



N 



which reproduces the exact result for the partition func- 
tion of N noninteracting spinless fermions. 

We now evaluate the contributions of the self-energy 
and of the gauge potential ao{i) up to the first non- 
vanishing order. This can be done in the usual way by 
making a loop expansion in the trace^Si^. Here we are in- 
terested in the lowest order contribution that survives in 
the low-energy and long- wavelength limit. This limit con- 
sists in expanding the one-loop contribution up to first 
order in dr and up to the second order in Ri — Rj im- 
plying that, eventually, we will set i j. This amounts 
to the so-called gradient expansion corresponding to the 
low-energy and long-wavelength limit of the action. We 
obtain in this way 

rrlog(l-G(o)za + G(o)S) = 

-Tr(G(o)*a) - Tr[i(G(o)SG(o)E) + ©(E^)]. (40) 

Note that rr(G(o)mG(o)S) = since T,^ = 0. Note that 
Ea. (j40p is invariant under global rotations. This imme- 
diately follows from the transformation law of the given 
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by Eq. (|27p accompanied by the similar transformations 



from the Ey 's: 



The a-dependent term in Ea. (|40|) contributes to the 
action in the following way 

Jo 

where G(o)i(0") lime^o G(o)i(-e), e > and 

Gm{r) = ^^-Oir)e-'^^. (41) 



The explicit representation (j4T|) tells us that 

Tr(G'(o)^a) ^ ©(e"^^), /i/3 » 1. (42) 
Let us now turn to the second term in Eq. (|40l) . We get 

-iTr(G(o)SG(o)S) 
= — ^ J G(Q)i{ti—t2)T.ij{t2)G(i^)j{t2~ti)Tiji{ti)dtidt2- 



2 ' ' 2 ' 



Introducing new variables, r 

and expanding the product ^ijiri + T)Y,ji{ri — r) = 
Sy(77)Sji(77) + 0(t) (this corresponds to the gradient 
expansion in imaginary time-^), this, to lowest order, re- 
duces to 

^\Y.f_ G(o){r)G^„^{~T)dTj\.,M^,,{r,)dv (43) 



With the help of Eq. (j4T]) we get 



-fi 



G(o)i'r)G^o){-T)d7 



a 1 

4 cosh^(/3^/2)' 



(44) 



Note that the energy scale is set by the chemical potential 
/i, so that the low-energy limit takes the form ^ 1. 

The effective spin action is then given by the sum of 
all the terms evaluated above: 



(45) 



where the SU(2) invariant measure is 

r, ( -\ TT dz^{t)dz,{t) 



and the effective action: 



. "'0 



f3 jeff 



(46) 



with the long- wavelength limit (i j) being implicit 
throughout our calculation. This action describes the 
SU(2) invariant ferromagnetic Heisenberg model with the 
effective coupling (see Appendix) 



eff 



Since 



J' 



eff 



ltyf)/(2cosh2^)<0. 



0(e-^^), fip » 1 



(47) 



this result tells us that the lowest order contribution to 
the free energy in E shows no magnetic ordering, pro- 
vided fj,(3 3> 1. However, one can safely truncate the 
expansion of the effective action (40) at second order, 
provided \t\f3 <^ 1. This, however, does not generate a 
high-temperature expansion of the free energy. The point 
being that the temperature is supposed to be " high" com- 
pared to the overall energy scale factor t. However, it still 
might be low compared with the "intrinsic" energy scale 
which is set by /i. Physically, the limit ^ 1 corre- 
sponds to a very small hole concentration^!. Therefore, 
our result (|T7)) provides us with a limited piece of infor- 
mation concerning the thermodynamic instability of the 
Nagaoka phase at finite temperature. It merely indicates 
that the FM order exponentially decays away from half 
filling at any finite temperature T = ^ 

To address the issue of the thermodynamic instability 
of the Nagaoka phase at any T > one should go back 
to the full series in the low-energy long-wavelength ex- 
pansion of the fermionic determinant (j40p and analyse 
further the dynamics of the spin variables as well. On 
the bipartite lattice with the nn interaction, that series 
can be summed up to yield the following contribution to 
the low-energy effective spin action: 



AS 



2/3 



Tr cosh(2/3E(i)), » 1. (48) 



The next step is now to calculate the long-wavelength 
asymptotic of this representation, which is not a trivial 
follow up and it still is in progress. Note only that this 
result is strongly dependent on the space dimensionality. 
However, there exists a universal feature of the U = oo 
Hubbard model phase diagram that naturally emerges 
from Eq. (|48p : the paramagnetic state is thermodynami- 
cally stable at any finite temperature at ^ = 0. This can 
be derived by considering the limit /i ^ oo in Eq. (|48|) 
at finite values of (3 and t. This result agrees with the 
qualitative arguments presented inS2,. In contrast with 
that, the MP treatment continue to predict in this case 
a fully polarized PM state. 



V. CONCLUSION 

In this paper, we formulate the improved SP MP the- 
ory to describe the PM phase of the Nagaoka limit of 
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the Hubbard model. It is clear that the physics behind 
the Nagaoka phase is controlled by strong electron corre- 
lations. Those correlations are in turn encoded into the 
NDO constraint. We improve the standard MF approach 
by taking the NDO constraint rigorously into considera- 
tion prior to the MF approximation. Once this is done 
we integrate out the fermionic degrees of freedom under 
the condition that the spin degrees of freedom are con- 
sidered at the MF level. In this way we arrive at the 
one-loop MF theory of the Nagaoka phase that predicts 
a FM ordering over the whole doping range. This result 
agrees with earlier variational estimates which take the 
NDO constraint fully into account but leaves some of the 
spin degrees of freedom completely frozen. At the same 
time, the conventional SF MF theory, that treats the 
constraint at the MF level, predicts the FM phase over 
a finite though quite large doping range. This result of 
the standard MF theory thus appears just as an artifact 
of the uncontrolled treatment of the NDO constraint. 

We show that the SF MF theory automatically drives 
the system into the FM ordered phase, and hence it can- 
not describe the underlying physics even in a qualita- 
tively correct way. This happens because of the fact that 
the spin degrees of freedom are only considered quasi- 
classically. However, it is the quantum spin dynamics 
that ultimately determines the magnetic properties of the 
Nagaoka phase. Therefore, to get any physically reason- 
able result one should enforce both the NDO constraint 
rigorously and treat the dynamics of the spin degrees of 
freedom beyond the MF approximation. 

As a first step in this direction, we derived the trun- 
cated low-energy long-wavelength effective spin action 
with the spin variables being treated beyond the MF ap- 
proximation. This approximation corresponds to the ex- 
pansion of the effective free energy in powers of a small 
parameter \t\f3 close to half filling. It turns out that treat- 
ing the spin dynamics beyond the MF approximation 
completely destroys the FM order predicted by the SF 
MF theory. Despite of that, this result does not prove 
the thermodynamic instability of the Nagaoka phase. To 
address this issue one should go back to the full non- 
truncated low-energy long-wavelength effective spin ac- 
tion as discussed in the previous sections of the paper. 
Nevertheless, we can add with certainty that the SF 
MF approach produces spurious results and is therefore 
not reliable for the description of the Nagaoka ferromag- 
netism. 



Acting with the "lowering" spin operator S~ on the 
"highest weight" state | t) we get the normalized su(2) 
CS parametrized by a complex number z 



^ATR 



:exp(z5-)| T> 



:(l T>+^U». 

(50) 



In the basis spanned by the vectors | t)j I i) we have 

5+ = |T)iai,^- = U)|(TI,^. = i(| T)I(T|-U)ia I). 

The CS symbols of the su(2) generators are then easily 
evaluated to be {S"' (z\S\z)): 



S' 



11 -Ul 



2 1 + |z| 



l + |z|2' 
^ 1/4, (S^U = 3/4. (51) 



There is a one-to-one correspondence between the 
su(2) generators (|49| and their CS (classical) symbols 
given by Eqs. (|5ip . Given a quantum Hamiltonian 
H — H{S), the corresponding imaginary time phase- 
space action takes on the form, 

d 

Asui2) (z, 2) = - (-^1 ^ + (52) 



with the kinetic term being given by 



1 zz — zz 



iz — z = -- 



'dt' ' 2 1 + Iz 
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In particular, for the quantum s = 1/2 Heisenberg model. 



one gets 



H^JJ2iS^S,~l/A), 
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VII. APPENDIX: SU(2) ALGEBRA AND 
COHERENT STATES 

Consider the su(2) algebra in the lowest s = 1/2 rep- 
resentation: 

[S,,S±]^±S±, [5+,5_]=25„ 5^2 = 3/4. (49) 



^From the geometrical viewpoint, the su(2) coherent 
states \z) can be thought of as sections of the magnetic 
monopole bundle P{S^, U{1)), with the U(l) connection 
one-form, ia, frequently refereed to as the Berry connec- 
tion. Base space of that bundle, two-sphere 5^, appears 
as a classical phase-space of spin, whereas its covariantly 
constant sections, \z) : 'Dt\z) :— {dt + ia)\z) — 0, form a 
Hilbert space of a quantum spin. 
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